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Abstract 

Some results about phase separation in coupled map lattices satisfying a conserva- 
tion law are presented. It is shown that this constraint is the origin of interesting 
antiferromagnetic effective couplings and allows transitions to antiferromagnetic and 
superantiferromagnetic phases. Similarities and differences between this models and 
statistical spin models are pointed out. 
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1 Introduction 



Coupled map lattices (CMLs) are spatially extended dynamical systems that 
have been considerably investigated recently on various grounds [1] . They are 
approximations to continuous systems particularly well suited for numerical 
and analytical calculations. From another point of view they can be considered 
as phenomenological models describing the dynamics of a large number of 
macroscopic chaotic subsystems. In both cases one is interested to study their 
behaviour at length scales larger than the scale where chaos is present. 

Particular attention has been dedicated to the study of collective behaviour in 
CMLs in presence of conservation laws. This class of models are considered rel- 
evant to describe several physical phenomena like surface waves in a container 
or disturbances in fluids, where mass, momentum and energy are conserved. 
A motivation for their interest is that, as it can be easily argued, the existence 
of conservative quantities should play an important role in determining the 
long distance properties in such systems. Several authors [2-4] studied coupled 
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map models with conservation laws as possible models for chaos in extended 
systems, examining phase transitions from ordered to chaotic states that are 
driven by the map parameters or by the value of the conserved quantity. 

CMLs have been studied [5-8] also from another point of view. In this case au- 
thors try to transfer concepts and results from equilibrium statistical mechan- 
ics to systems, like CMLs, whose dynamics is microscopically irreversible and 
does not satisfy detailed balance. In particular, studying systems of chaotic 
maps that exhibit a Ising - like symmetry, it can be shown that they undergo, 
in particular conditions, a phase - ordering dynamics. At large scale lengths 
the evolution of these deterministic systems has coarsening properties similar 
to Ising models or their continuous versions, the time - dependent Ginzburg - 
Landau equations [9,10]; the local chaos plays the role of the stochastic noise 
present in a heat bath or in the noise term of the Langevin equation. One 
expects obviously that, in the conserved order parameter case, this kind of 
CMLs should show a behaviour similar to the so called model B. This model 
is described by the Cahn - Hilliard [11] equation, leading to asymptotic state 
in which the phases occupy two large domains separated by a single boundary. 
Phase separation in CMLs with conserved dynamics was observed [12] in the 
case of a lattice of chaotic maps in contact with a thermal bath evolving with 
Kawasaki dynamics. In this case the temperature plays the role of an addictive 
noise, while chaos produces chaotically fluctuating couplings. It was deflnitely 
shown by J. Kockelkoren and H. Chate [13] that the coarsening process for 
this model is strictly similar to that of the corresponding Ising model. In the 
same article the authors proposed to study conserved dynamics of CMLs using 
the approach introduced by Y. Oono and S. Puri [14]. In this approach the 
dynamics corresponds to the discretization of the Cahn - Hilliard equation, 
each map on a site of the lattice representing the effect of a coarse - grained 
free energy. It was also shown that this model undergoes a phase transition 
between a short range and a long range ordered phase. 

The purpose of this paper is to extend the analysis of this model and to study 
the effect on the phase separation process of the coordination number used to 
discretize Laplacian operators on the lattice. In particular, it will be shown 
that the effect of the second Laplacian derived from the Cahn - Hilliard equa- 
tion, assuring the order parameter conservation, amounts to the introduction 
of an effective antiferromagnetic coupling. The competition between ferromag- 
netic and antiferromagnetic couplings in these systems plays an important role 
in the phase separation process. 
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2 The model 



Following [13,14] we consider a two-dimensional square lattice of coupled iden- 
tical maps / acting on real variables Xi. The discrete - time dynamics is gov- 
erned by the following equations: 



xr=^(x*)--^E^(^(4)-4), (2) 

where Af is the chosen number of neighbors of site i, the sum is over these 
neighbors and g is the coupling strength; periodic boundary conditions are 
assumed. The conservation of the order parameter 

M^J2^i (3) 



is incorporated by eq. (2), representing the second Laplacian in the Cahn 
Hilliard equation. 

The map used in the numerical simulations is the following: 



-| exp [a{x + ^)] if X E [—oo, — |], 

fix if X e [-|, |], 

^ exp [q;(| — x)] if xe[\,+oo\, 



(4) 



i.e. a modified version of the map used in [13], that was defined on the interval 

[—1,1]. The map 4 is defined for every x E R; the modification is motivated 
by the fact (ah^cady stressed in [4] and verified by the author) that, due to the 
redistribution step of the Oono - Puri dynamics (2), variables Xi{t) are not 
constrained to take value in [—1, 1]. Details on this map can be found in [15], 
where it was used for similar motivations. Choosing /j, — 1.9 and a — 6, f has 
two symmetric chaotic attractors, one with a; > and the other with a; < 0; 
this allows the unambiguous definition of Ising spin variables cr* = sgn[x*] 
associated to each dynamical system. 

To study the phase separation process, uncorrelated initial conditions were 
generated as follows: one half of the sites were chosen at random and the 
corresponding values of x were assigned according to the invariant distribution 
of the chaotic attractor with x > 0, while to the other sites were similarly 
assigned values with x < 0. With a good approximation the order parameter 
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M vanishes. We associated an Ising spin configuration {si{t)} = {sgn[xi(t)]} 
to each configuration of the x variable. Lattices from 256 x 256 up to 512 x 512 
with periodic boundary conditions were used. The average domain size R{t) 
was measured by the relation C[R{t),t] = 1/2, where C{r,t) = 
is the two point correlation function of the spin variables. R{t) was averaged 
over many different samples of initial conditions. More complicate correlation 
functions and related lengths, that will be introduced in the following, have 
been measured. Another variable is often considered in growth processes: the 
persistence p{t) [16], defined as the fraction of sites that have not changed their 
initial s values. It has not been evaluated because, to get a reliable calculation 
of persistence, one should use very large lattices which is not compatible with 
a fine scan of the coupling variable. 

One normally expects that, for large couplings, the conservation of the total 
"Magnetization" (2) leads to equihbrium states where there are only two large 
domains with aligned spins. However this conservation law is also compatible 
with more complicate phases. In statistical spin models, they are generated by 
the presence of repulsive couplings in the hamiltonian. In the sequel it will be 
shown that eq. (2), in addition to the conservation law, generates these new 
couplings. 



3 Nearest neighbors simulations 

As a first step we considered maps interacting with their nearest neighbors, 
corresponding to A/" = 4 in (1, 2). For various values of g the characteristic 
length R was measured as a function of time: R saturates for weak couplings 
at values small compared to the lattice size. For larger couplings it shows scal- 
ing behaviour and one gets complete phase separation in the spin variables. 
The value of g that discriminates between these two regimes is (7 ~ 0.03; 
its precise evaluation is beyond the purposes of this paper and will be pre- 
sented elsewhere. Figures 1 and 2 illustrate these changes. Fitting late times 
growth in the phase separation region with the law R{t) ~ At^ one gets the 
value z — 0.34 ± 0.01, in agreement with the result found by [13] in the 
case of 8 neighbors and with the Model B class of universality expectation 
(Lifshitz-Slyozov law). For higher values of the coupling one could expect a 
faster transition to the coarsening regime; however one finds a metastability 
region between g = 0.07 and g = 0.125. Studying i? as a function of time one 
finds at small values of a plateau whose extension grows with g. Between 
g — 0.120 and g — 0.125 the system gets trapped into blocked configurations 
with interfaces between equal phase domains completely pinned. At g = 0.127 
the behaviour of R(t) suddenly changes, the plateau being reduced to a fiex 
point; eventually, for 51 > 0.128 in few time - steps the lattice reaches a com- 
pletely ordered antiferromagnetic state (see fig. 3(a)). Antiferromagnetic (AF) 
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Fig. 1. Snapshots of a 256 x 256 sites CML with four neighbors for g = 0.01 (a) and 
g = 0.05 (b) after 10^ time steps. The full range of maps' values has been colored 
with a 16 gray levels scale from black to white. 
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Fig. 2. Time evolution of the domain size R{t) for g = 0.01 (a) and g = 0.05 (b). 
The solid lines is a best fit to the late time growth with power law R{t) = At^ with 
z = 0.341. 

ordering at large 51 is a normal phenomenon: for > a map at time t + 1 

gets a contribution with changed sign from its value at time t and configu- 
rations with aligned spins become unstable. This happens also in the case of 
CMLs with non conserved order parameter, i. e. in absence of the second step 
of the dynamics. However, in the present case the AF ordering is extremely 
precocious. The reason for this behaviour is the second redistributive step (2) 
of the conservative dynamics. Let us focus our attention on site i and let us 
suppose, for example, that the first step (1) increases, in average, the values 
of its neighbor maps. It is evident that, in the second step (2), this increment 
will result in a negative contribution to the Xi variable. This fact produces 
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Fig. 3. The antiferromagnetic ordering ai g = 0.128 for a CML with four neighbors: 
a) a magnified portion (100 x 100) of a 512 x 512 lattice after 2000 time steps, b) 
the probabihty distribution function for the same lattice. The full range of maps' 
values has been colored with a 16 gray levels scale from black to white. 

an effective antiferromagnetic coupling. As we said, at g = 0.128 the system 
quickly reaches antiferromagnetic ordering, and this phenomenon is preceded 
by a blocked phase and by met ast ability, which can be considered as pre- 
transitional effects. In fact, for lower values of g, one can see that the domains 
of aligned spins have a checkerboard structure in which values of maps be- 
longing to various peaks of the asymptotic probability distribution functions 
(PDFs) alternate (see fig. 4). The transition to the AF ordering occurs be- 
cause, for g > 0.128, the asymptotic PDF reduces to two narrow peaks of 
opposite sign (see fig. 3(b)). The AF ordering appears suddenly, the antiferro- 
magnetic correlation function displaying a jump after few tens of time steps. 
This ordering process survives for higher values of g; subsequently the dy- 
namics slows down until, by increasing g, the system evolves towards blocked 
configurations. As many characteristics of this behaviour are similar to the 
case of jV' = 8, we defer a more detailed analysis and discussion. Indeed, a 
new interesting phenomenology comes out from the study of this same system 
when "laplacians" are discretized using a higher number of neighbors, and it 
confirms the role of the second step of the dynamics as generator of an effective 
AF coupling . 



4 Nearest and next to nearest neighbors simulations 

For small g the behaviour of the CML in the case of eight neighbors is not 
substantially different from the case of four neighbors, apart from the fact 
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Fig. 4. Pre-transitional checkerboard structure ai g = 0.127: a) a magnified portion 
(100 X 100) of a 256 x 256 lattice after 10^ time steps, b) the probability distribution 
function for the same lattice. The full range of maps' values has been colored with 
a 16 gray levels scale from black to white. 

that the scale of g is smaller. For g < 0.007 no phase separation is attained 
and the system gets blocked when equal spin domains reach the dimension 
of few lattice spacings. For g between 0.007 and 0.110 the growth of opposite 
phases domains does not stop and one measures a growth exponent of | at 
late times. However, for intermediate times scales, in the interior of equal spin 
domains one can note the formation of a new sub arrangement of the maps, 
which corresponds to stripe-ordering. The stripes in each domain correspond 
to peaks in PDFs belonging to the same attractor of the map (4). Fig. 5 
shows for g = 0.065 a snapshot of the lattice and the sites PDF after 100000 
time steps. For g < 0.110 the competition between ferromagnetic and lamellar 
ordering processes is resolved in favor of the former, but this situation reverses 
for higher values of g. In this case the dynamics is essentially composed by two 
stages: starting from a random configuration (with zero average) of the map 
variables there is initially the formation of little striped domains oriented in 
the two directions of the lattice, then the motion of the domains walls produces 
a coarsening process. 

To study the growth of these striped domains we have considered, following 
[17,18], a new order parameter and two related correlation functions. These 
articles study the growth processes of striped domains in the superantiferro- 
magnetic (SAF) phase of the Ising model with nearest and next to nearest 
neighbors interactions. We recall that the SAF phase, in which the ground 
state is four-fold degenerate and consists in alternate up and down spins rows 
or columns, is related to the existence of an antiferromagnetic coupling be- 
tween next to nearest neighbors. In fact, if we call Ji and J2 the couplings 
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Fig. 5. The underlying striped structure for a CML with eight neighbors g = 0.065: 
a) a magnified portion (100 x 100) of a 256 x 256 lattice after 10^ time steps, b) 
the probability distribution function for the same lattice. The full range of maps' 
values has been colored with a 16 gray levels scale from black to white. 

between nearest neighbors and between next to nearest neighbors respectively, 
the SAF phase corresponds to | Ji| < 2| J2I and J2 > 0. The existence of striped 
domains is therefore another corroboration for the existence of an effective an- 
tiferromagnetic coupling in the model here presented. Let us divide the lattice 
in 2 X 2 cells and consider a two-components local order parameter 

^'^^ ^'1 (5) 



defined in each cell a in the following way: 

„i,S _.a I _.S _.a _.S „l3 _.a _q _.q i _.S //^\ 



where erf are the clockwise ordered spins of the cell a. This order parameter 
allows a univocal labelling of the four SAF phase ground states. * allows the 
introduction of two correlation functions [18]: 



r£(r, t) = -iijfijt^^) + -i^P^^Pt"^'), (7) 
rt(r, t) = lii^fi^f^'-') + lii^^i^t'''), (8) 

where x and y are unit vectors in the x and y directions and (...) indicates the 
average over lattice cells and different initial conditions. and Ft measure 
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Fig. 6. a) The three measured lengths R (circles), Ri (squares), Rt (triangles) as a 
function of time at g = 0.105; the solid lines is a best fit to the late time growth 
with power law R{t) = At^ with z = 0.334. b) The same quantities at 5 = 0.111; a 
power law fit gives in this case z = 0.337 for Ri and z = 0.326 for Rt. 

respectively the correlation properties in the direction where the spins are 
aligned and in the direction where the spins are alternate. 

Similarly to C{r,t), it is possible to get from the measurement of ri{r,t) 
and rt{r,t) two characteristic lengths Re{t) and Rt{t) and to point out the 
competition between the ordering processes. Figure 6(a) shows R, Re and Rt 
as a function of time Bit g — 0.105. There is an initial increase of all these 
three lengths, then, after few dozens of time steps, the ferromagnetic ordering 
prevails. Re and Rt go to zero and the growth of R goes on, reaching the 
scaling regime power law with exponent |. For higher values of g the situation 
reverses. Fig. 6(b) shows what happens at g — 0.111; Re and Rt get quickly 
the Lifshitz-Slyozov law regime. These two evolutions are also represented in 
fig. 7 and fig. 8 through two snapshots of the lattice at t = 100 and t — 50000 
for these two values of g. Another interesting feature is the presence of an 
anisotropic growth of domains in the SAF phase similar to that measured 
by [18]: Re is always greater than Rt, as long as we are far enough from 
the formation of a single domain covering the whole lattice. This happens, as 
reported in [18], when the coupling between nearest neighbors is ferromagnetic, 
while the reverse is true when this coupling changes sign. 

This situation, in which one has complete striped ordering, persists increasing 
g until gi ~ 0.13. For higher values the incapacity of the system to get complete 

phase separation couples to a behaviour of the maps' asymptotic PDFs that 
have characteristics similar to the case A/" = 4. In fact, while at ~ 0.13, e. g., 
PDFs display two narrow peaks at opposite values of the x variable (fig. 9(a)), 
increasing g they display first the broadening of these peaks, then the birth 
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Fig. 7. Snapshots of a 256 x 256 sites CML with eight neighbors for g = 0.105 (a) 
after 100 and (b) after 50000 time steps. The full range of maps' values has been 
colored with a 16 gray levels scale from black to white. 




Fig. 8. Snapshots of a 256 x 256 sites CML with eight neighbors for g = 0.111 (a) 
after 100 and (b) after 50000 time steps. The full range of maps' values has been 
colored with a 16 gray levels scale from black to white. 

of two new peaks at symmetric x positions (fig. 9(b)). A further increase of g 
gives rise to an increasing number of peaks (fig. 9(c)). 

At this point, order parameters like those borrowed from the Ising model, 
the ferromagnetic, antiferromagnetic and lamellar ones, become inadequate 
to describe the dynamics of these systems. Even when the lengths related 
to these order parameters have a monotonous behaviour, one finds growth 
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Fig. 9. Probability distribution functions for CMLs with eight neighbors after 10^ 
time steps for g = 0.13 (a), 5 = 0.20 (b), 3 = 0.30 (a). 

exponents different from the expected value of |. This exponent is related 
to the mechanism decreasing the curvature between domains in which the 
order parameter has an homogeneous value and to the conservation law. It is 
worthy to mention that the systems we are studying satisfy, in this regime, 
the conservation law of a quantity that is no longer the true order parameter 
for the coarsening process. 



5 Conclusions 



In this paper some results on the time evolution of lattices of coupled maps 
(to which an Ising spin can be associated) in presence of a conservation law 
have been presented. Studying both the cases of nearest and next to nearest 
neighbors, intervals of the coupling have been located in which the coarsening 
process is similar to the classical growth phenomena described by model B, 
confirming that also deterministic systems like CMLs are able to show similar 
processes of phase separation. It has been pointed out an interesting feature, 
i.e. the effective antiferromagnetic coupling rising from the conservation law 
constraint. Normally the study of these systems has been limited to low values 
of the coupling g and antiferromagnetic effects have not been considered. The 
present analysis shows that they give rise to an interesting phenomenology, 
including a scaling region and a phase diagram which resembles in some case 
the Ising model. For example, for A/" = 8, starting from a ferromagnetic phase, 
one can go, increasing the coupling value, to a Superantiferromagnetic phase 
through a paramagnetic one. However, particularly for strong couplings, we 
have shown the greater complexity of these models. This complexity cannot 
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be described by concepts and quantities transferred from the study of spin 
statistical models. 

Furthermore this study confirms the importance of theoretical investigations 
of the asymptotic probability distribution functions [19,20] with respect to 
the ordering process. With regard to this, one cannot exclude that, after a 
time longer than the one used in these simulations, eight neighbors CMLs 
evolve definitely towards SAF configurations and their PDFs towards two 
peaks structure also in the case of strong coupling. 

Finally, in this analysis the parameter in (2), corresponding to the mobility 
coefficient in the Cahn - Hilliard model, has the fixed value l/A/", so that, with 
the subsequent sum, makes an average of the neighbors' value increments. Ex- 
ploring the role of this parameter could be the subject of future investigations 
on these systems. 
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